机 械 工程 学 dk 
JOURNAL OF MECHANICAL ENGINEERING 


DOI: 10.3901/JME.20* 5 443 


F ET ZEE LAM. HJ ER [8] e SARH RIETTE 
成 算法 
谢 展 阳 1 Xe? 肖 光 明 3 张 贤 杰 ! IAE! FIF 
(1. 西北 工业 大 学 机 电学 院 ”西安 710072; 


2. 中 国航 空 工业 集团 公司 雷 华 电子 技术 研究 所 ”无锡 214063; 
3. 中 航 西 安 飞 机 工业 集团 股份 有 限 公 司 西安 710089) 


摘要 : 提出 一 种 单身 纤维 增强 复合 材料 (UD-FRP) 周 期 性 蜗 保 真 代表 性 体积 元 (RVE) 的 创建 算法 。 首 先 基于 随机 游 走 算法 ， 采 
用 贝 塞 尔 (Bézier) 曲 线 表 征 纤 维 形态 , 使 曲线 各 控制 点 处 的 切线 方 回 的 统计 分 布 服从 多 变量 冯 米 斯 - 引 租 尔 (Von-Mises Fischer) 
函数 ， 创 建 出 初步 的 软 核 模型 。 进 而 基于 力 偏 置 算法 ， 在 保证 RVE 内 部 纤维 不 产生 过 度 扭 曲 的 前 提 下 ,调整 纤维 间 的 重 半 ， 
创建 出 对 应 的 硬 核 模 型 。 调 整 过 程 中 ， 通 过 引入 三 维 近邻 列表 算法 ， 大 幅 提高 算法 效率 。 最 后 对 所 创建 的 盒 有 纤维 波纹 度 特 
征 的 RVE 模型 进行 了 横向 剪 切 渐进 损伤 分 析 ， 用 以 揭示 该 方法 的 研究 价值 及 应 用 领域 。 结 果 表 明 ， 该 方法 能 够 快速 有 效 地 
创建 不 同 尺 寸 及 纤维 波纹 度 的 RVE 几何 构 型 , 这 种 含有 波纹 度 特征 的 纤维 构 型 会 对 UD-FRP KIPR I8] BIH PE pEr E — XE 86 Un, 
纤维 波纹 度 越 大 ，UD-FRP 的 横向 剪 切 强度 越 大 。 
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Abstract: A generation algorithm of periodic high fidelity representative volume element (RVE) model of unidirectional 
fiber-reinforced composites (UD-FRP) is introduced. Firstly, based on the random walk algorithm, the fiber morphology is 
characterized as the Bézier curve, and the tangent direction at each control point is adjusted to follow the multivariable von Mises 
Fischer function to generate a preliminary soft-core fiber network. Then based on the force-biased algorithm, the corresponding 
hard-core model is created by adjusting the overlap between fibers to ensure that the fibers in RVE do not twist excessively. In the 
adjustment process, a three-dimensional near-neighbor list is introduced to significantly improve the algorithm's speed. Finally, the 
RVE containing fiber waviness is preliminarily calculated for transverse shear, so as to reveal the research value and application field 
of this method. The results show that the RVE of the required size and fiber waviness can be generated effectively and quickly by this 
method, the transverse shear strength of UD-FRP can be affected by the fiber configuration with waviness characteristics. The 


transverse shear strength of UD-FRP increases with the increase of fiber waviness. 
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0 前言 


纤维 增强 聚合 物 复合 材料 FRP) 特殊 的 比 刚 
度 和 强度 使 其 成 为 结构 应 用 的 理想 材料 ， 其 结构 件 
在 航空 航天 、 海 洋 、 土 木工 程 、 汽 车 等 行业 的 应 用 
日 益 广 泛 。 在 复合 材料 制造 过 程 中 ， 预 浸 料 或 干 布 
中 的 纤维 构 型 通常 并 不 是 理想 的 直 纤 维 ， 对 于 单 根 
纤维 而 言 ， 纤 维 的 局 部 错位 导致 该 处 主 方 回 与 整 根 
纤维 的 预期 主 方向 存在 一 定 偏差 ， 这 种 现象 通常 被 
称 为 纤维 波纹 度 山 。 这 种 纤维 本 身 构 型 特点 是 导致 
固化 过 程 中 产生 残余 应 力 的 原因 之 一 ， 同 时 也 会 对 
成 型 后 的 复 材 件 的 力学 性 能 造成 显著 影响 P31。 

纤维 的 波纹 度 特 征 往往 会 导致 纵向 压缩 载荷 下 
的 纤维 微 届 曲 由 ， 这 正 是 纤维 增强 复合 材料 主要 损 
伤 机 制 之 一 ， 并 且 这 种 微 届 曲 也 使 得 复合 材料 的 弯 
曲 及 前 切 等 力学 性 能 产生 较 大 波动 ， 这 种 波动 则 会 
导致 复合 材料 件 力学 性 能 上 的 不 确定 性 名。 因此 ， 
准确 表征 纤维 的 这 种 波纹 度 特征 具有 重要 意义 。 

为 准确 表征 FRP 失效 过 程 中 的 力学 行为 , 有 必 
要 建立 微观 结构 的 高 保 真 代表 性 体积 元 CRVEO H 
于 分 析 其 均 质 化 特性 。 目 前 已 有 大 量 的 工作 致力 于 
单 向 纤维 增强 复合 材料 CUD-FRPO 高 保 真 微观 模 
型 的 构建 ， 但 其 中 大 多 数 仅 涉及 二 维 RVE 几何 模 
型 , 聚焦 于 定义 RVE 内 部 的 各 种 纤维 分 布 特征 和 生 
成 算法 研究 ， 如 探寻 法 及 随机 顺序 吸附 CRSAO 算 
法 67。 三 维 的 RVE 模型 则 是 通过 沿 纤 维 方 同 简 单 
地 拉 伸 二 维 RVE 得 到 ， 所 有 纤维 均 是 平 直 排 列 。 

现 有 研究 表明 ， 纤 维 的 波纹 度 特征 可 以 通过 解 
析 式 表达 (规则 波纹 度 ), 也 可 通过 随机 过 程 来 确定 
纤维 的 几何 特性 (随机 波纹 上 度 )。 目前 很 省 有 天 然 纤 
维和 普通 纤维 实现 随机 波纹 度 特征 的 研究 报 
告 ， 若 创建 出 可 调控 几何 结构 特征 (纤维 分 布 和 波 
纹 度 ) 的 RVE 模型 , 便 可 进一步 用 于 实际 微观 结构 
的 均 质 化 分 析 。 

在 微观 力学 计算 中 ， 数 值 均匀 化 需要 使 用 适当 
的 周期 边界 条 件 (PBC)，RVE 边界 条 件 施加 的 合 
理性 直接 决定 了 模型 的 计算 效率 和 精度 ， 进 而 影响 
模型 的 均匀 化 参数 和 相关 机 械 性 能 的 计算 准确 度 
UU, qf PBC 的 前 提 即 要 求 RVE 具有 几何 的 周 
期 性 。 


为 了 描述 纤维 波纹 度 特征 对 UD-FRP 轴 问 损伤 
特性 的 相关 影响 ，Catalanotti 和 Sebaey 等 人 开发 了 
一 种 以 统计 方式 表征 纤维 波纹 度 的 半 随 机 算法 
L213]。 该 算法 首先 将 连续 纤维 离散 为 多 组 球体 ， 通 
过 电子 计算 机 断层 扫描 (CT) 确定 UD-FRP 内 部 各 
组 纤维 球 的 伺 转 角度 ， 根 据 冯 - 米 塞 期 (von-Mises ) 
分 布 对 CT 数据 进行 统计 拟 合 ， 以 浓度 参数 表征 纤 
维 波 纹 度 的 大 小 ， 最 后 利用 随机 扰动 程序 及 优化 程 
序 保证 RVE 与 实际 试 件 具有 相同 的 浓度 参数 , 同时 
保证 了 模型 自身 的 几何 周期 性 。 该 算法 虽 构 建 出 与 
实际 UD-FRP 在 统计 意义 上 具有 相同 纤维 波纹 度 的 
RVE 模型 , 但 随机 扰动 程序 仅 单独 作用 于 各 纤维 球 
上 ， 没 有 考虑 纤维 上 相 邻 球体 间 的 相互 作用 ， 导 致 
纤维 在 局 部 可 能 出 现 过 度 的 扭曲 ， 这 在 实际 的 纤维 
结构 中 是 不 存在 的 。 

针对 上 述 研究 和 分 析 ， 可 见 构建 包含 波纹 度 的 
周期 性 RVE 可 以 准确 用 于 探究 波纹 度 特征 对 
UD-FRP 失效 行为 的 影响 ， 但 是 需要 考虑 模型 目 号 
的 物理 真实 性 。 本 文 针 对 这 一 问题 ， 提 出 了 一 种 基 
于 随机 游 走 算法 M4151, 考虑 RVE 几何 周期 性 及 纤维 
波纹 度 特征 的 UD-FRP 高 保 真 微观 力学 模型 生成 算 
法 ， 并 将 其 用 以 损伤 失效 的 相关 有 限 元 计算 ， 这 对 
解决 FRP 的 实际 应 用 问题 具有 重要 意义 。 


1 算法 架构 


考虑 纤维 波纹 度 的 UD-FRP 周期 性 RVE 创建 主 
要 由 两 个 步骤 构成 : 软 核 系 统 构建 及 便 核 系统 调整 ， 
如 图 1 所 示 。 

软 核 系统 构建 过 程 首 先 根据 RVE 纤维 半径 、 纤 
维 体积 分 数 及 RVE 尺寸 大 小 的 要 求 , 采用 随机 扰动 
方法 创建 出 纤维 的 二 维 随机 位 置 点 ， 其 创建 方法 已 
有 诸多 文献 报告 ， 此 处 不 再 著述 16171。 随 后 沿 纤维 
主 方 同 采用 随机 游 走 算法 在 已 有 的 纤维 位 置 点 依次 
构建 三 维 纤维 构 型 ， 并 最 终 对 整个 纤维 系统 进行 几 
何 周期 性 设置 。 该 过 程 允许 纤维 间 的 相互 重 铬 ， 因 
而 称 为 软 核 纤维 系统 。 

便 核 系统 调整 过 程 则 是 在 已 构建 的 软 核 纤 维系 
统 基础 上 进一步 加 以 调整 ， 在 纤维 构 型 基本 不 变 的 
前 提 下 ， 消 除 纤维 之 间 的 重 骆 现象 。 该 调整 过 程 引 
入 力 偏 置 算法 ， 即 定义 纤维 间 斥 力 / 引 力 - 相 对 距离 
及 纤维 自身 弯曲 力 -曲率 的 关系 式 , 并 将 其 应 用 于 已 
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创建 的 软 核 纤 维系 统 。 此 外 ， 为 加 速 调 整 效 率 ， 引 
入 三 维 近 邻 列 表 算 法 , 该 算法 由 William Mattson 等 
人 给 出 具体 实现 细 市 1181。 
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图 1 算法 框架 示意 图 


2 重 登 纤维 软 核 系 统 创建 


2.1 随机 游 走 算法 生成 单 根 纤维 

根据 RVE 纤维 半径 及 纤维 体积 分 数 的 要 求 , 可 
简单 创建 出 随机 的 三 维 纤维 位 置 点 ， 为 创建 包含 波 
纹 度 特征 的 纤维 构 型 ， 利 用 纤维 随机 游 走 算法 依次 
单独 创建 出 每 一 根 纤 维 , 该 过 程 中 采用 Bézier 曲线 
表征 每 根 纤 维 的 形态 ， 曲 线 各 控制 点 的 方 同 均 服从 
iU 3E Hr-9* x 7K«(Von-Mises Fischer) 分 布 。 访 分布 
作为 一 种 第 见 的 分 布 水 数 ， 和 用 以 描述 名 望 共 有 优 
先 方 回 以 及 相应 可 靠 性 参数 的 目标 方 辐 ， 如 下 式 所 
A 


kuu 
f| n. o I (1) 
XX u Z3 Hs Imi&. URREN, k 
则 为 相对 应 的 可 对 性 参数 。 
本 算法 中 各 根 纤 维 曲线 的 控制 点 个 数 y 由 RVE 
的 轴 加 长度 瑟 及 纤维 半径 六 共同 诀 定 : 
5 二 (17)+1， 并 定义 第 ii 二 1,.…,5) 个 控制 点 的 


7j |n] m 为 第 i 个 控制 点 与 第 1-1 个 控制 点 的 相对 方 
H, BẸ 7 = CN Yi T Yis Zi =R o 

Z5 ES SIE AH] ES RAE Ee UL ZK 3€ jp - 39 RC 
布 的 随机 游 走 算法 ， 将 导致 较 长 的 纤维 曲线 可 能 杰 
曲 成 圆 形 ， 不 符合 UD-FRP 的 几何 要 求 。 因 此 为 确 
保 所 创建 纤维 曲线 物理 结构 的 真实 性 ， 采 用 多 变量 
冯 米 堵 斯 一 费 舍 尔 分 布 以 保证 纤维 曲线 的 各 控制 点 
同时 与 UD-FRP £555] 77 [8] , BJ RVE 全 局 主 方向 和 
本 一 个 控制 点 的 切线 方 回 同时 具有 一 定 的 相关 性 。 
该 分 布 函 数 由 单 变 量 汉 米 起 斯 一 费 舍 尔 分 布 扩展 而 
来 ， 可 描述 同时 上 共有 两 个 优先 方 同 回 量 M, u, 以 及 
对 应 可 菲 性 参数 ,k, 的 目标 方 回 ， 如 下 式 所 示 : 
fu | B Ks uk) = Iu Rk, uk.) T cS (2) 


ku, +k, u 
I(u,, ku, k;)-— 一 一 一 


jM SMS) oy PuMIT SMS 
glet ku; | To +k ) 


上 式 中 第 一 个 优先 方 同 Wj =(0,0,1) 为 RVE 的 全 局 主 
方向 。 第 二 个 优先 方向 u, (x,y,z) 则 代表 每 一 个 控 
PRISAJE, BAXMA |a=; ki o 
别 为 其 对 应 的 相关 性 参数 。 对 比 可 得 ， 每 一 个 多 变 
量 VME 分 布 都 可 以 写成 市 有 参数 的 经 典 vMF 分 布 ， 
参数 由 下 式 确 定 : 


(3) 


k = ku, uA (4) 
ku, Tk 
u = uH : 2M3 (5) 


对 于 本 文 的 模型 ，ki 值 的 大 小 决定 纤维 各 控制 
点 方 同 与 RVE 全 局 主 方 回 的 相关 性 , k ERKA 
定 了 控制 点 与 前 一 个 控制 点 方 同 的 相关 性 ， 两 者 共 
同 决 定 纤 维和 目 遇 的 波纹 度 特 征 。 针 对 本 文 的 
UD-FRP, ?4 kı ~H, k 越 小 则 纤维 控制 点 游 走 
的 随机 性 越 高 ， 即 纤维 的 波纹 上 度 越 大 ， 刁 越 大 则 情 
况 相 反 。 


L 塞 尔 曲线 的 控制 点 


2 单 根 纤维 曲线 的 控制 点 空间 参数 
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在 描述 每 个 控制 点 的 具体 生成 过 程 之 前 ， 需 要 
明确 算法 中 各 控制 参数 , 以 RVE 中 的 第 7 根 纤维 为 
例 ， 如 图 2 所 示 。7 即 为 控制 点 的 间距 ， 由 上 文 可 
知 其 取 为 纤维 的 半径 。 yui(j 2L. MERRE 根 
纤维 的 轴 辣 方向 , 通常 取 为 RVE 的 全 局 主 方向 ， 访 
RVE 共 包 含 n 根 纤维 。 取 w(i=1,2,3) 为 控制 点 的 
方向 向 量 pu, 在 坐标 系 各 轴 上 的 投影 分 量 的 大 小 ， 
此 对 于 纤维 上 的 每 个 控制 点 的 方向 向 量 
Hy v Ww, w); F: 


w +w, +w =1 (6) 
对 于 每 根 纤维 曲线 ,所 示 算 法 的 具体 过 程 如 下 : 
Step 1: 为 保证 | 中 |=1， 且 服从 多 变量 vMF 分 

布 ,定义 : 


w -1—-M -N, (7) 
w,-N1-M-N, (8) 
w, =M (9) 


其 中 M 5 N 均 为 独立 的 随机 变量 ， 取 NN 为 单位 圆 
上 的 随机 数 ， 有 : 


N? +N;=1 (10) 
而 M 38 08) 4H ES SC: 
kv 
fec aD 
€ —e 
Step 2: 计算 vMF TAWAD ARAF): 
kx —k 
F(x) = L—— (12) 
e —e 


Step 3: 为 生成 符合 vMF 分 布 的 随机 变量 ， 计 算 累 
积分 布 函数 的 反 函 数 书 -(y) 如 下 式 : 

F's oge +A- 03) 
其 中 y 为 (0,1) 上 的 随机 数 ,代入 上 式 计算 之 后 使 得 
M =F (9); 
Step 4: 依次 计算 得 到 纤维 曲线 各 个 控制 点 的 方 问 
器 量 H5; 
Step 5: 给 定 RVE 模型 的 几何 尺寸 L(i=1,2,3)， 则 
每 根 纤 维 上 的 控制 点 个 数 得 以 确定 : 

s=(L,/r)+1 (14) 

Step 6: 计算 得 到 纤维 上 各 个 控制 点 的 坐标 ; 
Step 7: 生成 Bézier 曲线 ， 此 时 曲线 总 长 大 于 给 定 
RVE 的 高 上 度 L ， 曲 线 的 解析 式 如 下 : 


B(t) = X5 (1— £y ie (15) 


22 纤维 主 方向 与 球 链 坐标 标定 
上 文 算法 均 使 单 根 纤 维 生 成 于 原点 位 置 ， 并 且 


将 Bézier 曲线 沿 z 方 癌 平 均 赋 予 关 个 球体 ， 即 纤维 
由 一 串 连 续 球 链 构 成 ，r 的 值 只 需 保 证 最 终生 成 纤 
维 的 光 顺 性 即 可 ， 如 图 3 中 灰色 球 链 所 示 。 

虽然 初始 取 辐 u? =(0,0,1) 是 根据 全 局 主 方向 
进行 选择 的 , 但 由 于 该 算法 中 ji 的 随机 性 而 导致 最 
终 的 平均 纤维 取 癌 产生 偏 莽 。 为 了 保证 UD-FRP 的 
全 局 取 问 分 布 ， 旋 转生 成 的 的 纤维 曲线 ， 使 平均 纤 
维 取 问 与 RVE WEHEN IE JW =(0,0,1) 保 持 一 致 ， 


在 此 过 程 中 纤维 的 结构 (长 度 、 半 径 、 要 曲 ) 始 终 保 
持 不 变 。 
z (Zz) 


纤维 球 链 


图 3 单 根 纤 维 的 创建 示意 图 
e 为 1 Wu, 的 法 向 量 ，@ 为 J 与 JW 的 偏差 来 
fü. e 为 旋转 前 该 根 纤维 上 第 K(k = Lun) 个 球 
Jk r^ 相对 该 纤维 上 第 一 个 球体 x"” 的 方向 : 


æ = x" —x (i =1,2,3) (16) 
e= u,x ui (17) 
则 可 定义 旋转 矩阵 DD 为 : 
D(o,e,a) 
(18) 


= (e-0)e *- cosa((exmQ)xe)--sinac(e xo) 
其 中 “…” 表 示 标 量 积 ,“x? 表 示 又 积 。 因 而 纤维 上 第 
k 个 球体 旋转 后 的 坐标 为 其 旋转 定 阵 与 该 纤维 上 第 
一 个 球体 的 坐标 之 和 : 
x. =x -D(0,e,a) (19) 

最 终 将 旋转 后 的 纤维 依次 整体 平移 至 RVE 空 
间 内 的 相对 位 置 上 ， 通 过 该 随机 游 走 算法 对 每 一 根 
纤维 进行 上 述 操作 ， 即 可 初步 建立 了 一 个 由 重 郑 球 
链 构 成 的 软 核 纤维 模型 。 
2.3 几何 周期 性 设置 

在 复合 材料 多 尺度 力学 分 析 中 ， 使 用 周期 性 
RVE 可 大 幅 减 少 分 析 的 难度 及 计算 量 ， 且 数值 的 均 
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匀 化 通常 需要 施加 周期 性 边界 条 件 ， 为 满足 计算 所 
需 的 周期 性 边界 条 件 设置 ， 需 要 对 模型 划分 周期 性 
的 网 格 结构 , 前 提 则 是 所 建立 的 RVE 几何 模型 需要 
治 三 个 方 回 呈现 周期 性 的 特点 091。 

首先 对 所 建 纤 维 球 链 点 的 坐标 进行 位 置 判 定 ， 
上 文 已 将 每 根 代表 纤维 的 Bézier 曲线 使 用 一 组 等 间 
距 的 球体 进行 了 赋值 ， 并 且 通 过 旋转 矩阵 可 以 保证 
RVE WEHA C z 方向) 具有 周期 性 的 特点 ， 
因此 为 实现 另外 两 个 方 回 的 周期 性 , 现 将 RVE 的 所 
有 球体 按 z 方 同 以 球 链 间 距 工 (x -D 划 分 为 三 个 平 
面 PG=1,,r)， 单 独 对 每 一 个 平面 上 的 球 链 点 进行 
周期 性 条 件 的 设置 ， 如 图 4 所 示 ， 需 注意 图 4 中 仅 
以 纤维 体积 分 数 58.9%, RVE 三 个 方向 尺寸 均 为 
108 um 的 纤维 系统 作为 示意 概念 图 , RVE 具体 的 纤 
维 数量 则 由 初始 创建 的 二 维 纤维 坐标 点 数据 确定 。 

基于 周期 性 条 件 对 RVE 边界 处 的 几何 特征 坐 
标的 周期 性 特点 要 求 09， 对 于 各 平面 内 的 球体 ， 需 
要 考虑 n 个 球 心 ( 图 4 中 蓝 色 球体 ) 的 坐标 x (i=1,2; 
k=1,...,n) 超出 了 RVE 的 几何 范围 工时 ， 则 应 在 其 
对 称 位 置 生成 一 个 新 的 球体 以 保证 RVE 的 几何 周 
期 性 , 其 它 完全 位 于 RVE 范围 内 部 的 球体 则 无 需 改 


变 其 坐标 值 。 
Lu <r 
A =40,r <x < L -r (20) 
-L,x > L -r 
x 二 xX/ rA, (21) 


位 于 平面 四 个 角 点 处 的 球体 (图 4 中 红色 球体 ) 
应 对 应 其 他 三 个 角 点 处 的 球体 同时 保证 其 周期 性 ， 
即 在 另外 三 个 角 点 的 对 应 位 置 同 时 创建 对 应 周期 性 
条 件 的 球体 。 


内 部 球体 
MM ^ 边界 球体 
v 7^ 角 点 球体 


boo? 00906 9^ ~x 


图 4 RVE 周期 性 边界 条 件 设置 概念 图 
依次 对 地 个 平面 上 的 球体 进行 上 述 周 期 性 设置 ， 即 
可 保证 该 软 核 纤 维系 统 整 体 的 几何 周期 性 。 


3 非 重 登 纤 维 使 核 系 统 创 建 


3.1 三 维 近 邻 列 表 算 法 

上 文 已 创建 出 满足 所 需 纤 维 波纹 度 特 征 的 软 核 
纤维 系统 ， 该 系统 中 各 纤维 之 间 人 允许 发 生 重 登 ， 因 
此 在 非 重 于 纤维 硬 核 系 统 的 调整 过 程 中 需要 对 重 装 
的 纤维 球体 间 施 加 力 偏 置 算法 ， 以 分 离 相 互 重 车 的 
纤维 。 

为 施加 该 扰动 ， 准 确 判断 纤维 球 链 上 的 每 一 个 
球体 与 其 他 纤维 球 链 是 否 重 铸 及 重 车 的 球体 坐标 尤 
为 重要 ， 而 这 将 耗费 大 量 的 计算 时 间 。 本 文 引 入 三 
维 近 邻 列表 算法 ， 用 以 加 速 该 过 程 。 每 个 球 心 的 近 
邻 列表 中 包含 了 所 有 可 能 与 其 重 倒 的 其 他 球 心 点 
(距离 小 于 半径 之 和 ), 在 近邻 列表 中 搜索 目标 纤维 
球 链 上 所 有 球 心 的 近邻 球 心 ， 能 快速 找到 与 之 重 装 
的 纤维 球 链 以 及 产生 重 膨 的 小 球 。 在 施加 作用 力 之 
前 ， 可 以 快速 找到 可 能 与 之 重 闭 的 纤维 球 链 ， 然 后 
在 后 续 创 建 硬 核 纤 维系 统 的 过 程 中 只 需 考 虑 纤维 球 
链 是 否 与 其 近邻 列表 内 的 球体 充 登 ， 即 利用 三 维 近 
邻 列 表 的 方法 加 速 纤 维 球 链 间 重 又 情况 判断 的 过 
程 ， 提 高 本 算法 的 效率 。 

在 创建 近邻 列表 之 前 ， 将 空间 分 割 为 多 个 立方 
体 小 区 域 ， 如 图 5 所 示 。 小 立方 区 域 的 边 长 不 小 于 
纤维 半径 ， 其 每 个 顶点 均 设 为 结 点 。 区 域内 的 所 有 
球 心 点 的 坐标 同 下 取 整 后 分 配 到 对 应 的 结 点 处 
Cq eq. 便于 缩小 搜索 范围 。 取 每 一 个 球体 球 心 
点 为 其 三 维 近 邻 列 表 的 原点 (近邻 列表 的 中 心 点 )， 
创建 每 个 球体 的 三 维 近 邻 列 表 ， 则 近邻 列表 中 即 包 
括 所 有 可 能 与 原点 重 登 的 近邻 球体 ， 进 而 得 到 所 有 
与 该 球 心 可 能 重 登 的 结 点 。 通 过 仅 读 取 各 球体 近 令 
列表 内 的 其 他 球 心 点 ， 可 以 省 略 近 邻 列 表 外 球体 的 
影响 , 大幅 减 少 碰撞 判断 及 施加 作用 力 的 迭代 次 数 ， 
从 而 加 速算 法 。 


图 5 三 维 近 邻 列表 示意 图 


同时 可 将 决定 RVE 几何 周期 性 的 外 围 纤维 球 
链 和 内 部 其 他 纤维 球 链 的 近邻 列表 数据 分 开 存 储 ， 
自 完 对 外 围 纤 维 球 链 进行 便 核 系统 的 创建 过 程 ， 随 
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后 固定 其 球 心 坐标 ， 对 内 部 的 其 他 纤维 球 链 进行 该 
过 程 ， 以 此 保证 RVE 的 几何 周期 性 不 会 被 破坏 。 
3.2 纤维 则 斥 力 

由 于 要 求 RVE 中 的 纤维 与 纤维 之 间 不 可 产生 
重 登 ， 因 此 建 模 过 程 中 当 不 同 球 链 《〈 不 同 纤 维 ) 存 
在 重 登 时 ， 需 要 施加 斤 力 以 恢复 其 刚性 接触 。 需 要 
注意 的 是 ， 本 文 的 方法 中 人 允许 同一 球 链 内 部 的 球体 
BB. WK 6 所 示 ， 游 走 过 程 中 两 球 链 相 互 重 邯 ， 
其 半径 同 为 x 。 而 其 中 重 闭 度 最 大 的 两 球 ( 不 同 球 
EE 位 置 点 分 别 为 q 、9q,。 两 个 球体 间 的 斥 力 与 
其 重合 度 相关 , EA NA 为 两 球 互 相 重 只 的 区 域 大 


小 : 
Al 
图 6 两 球 重 登 模型 
Al - Pr — la, -q, | (22) 
则 球 q, 受到 的 斥 力 为 : 
Ep qd; 4; (23) 
© 2 la, -4;| 


AUD K, KEREI ASA, P ARAE CE 
ERKE HSJLUSR GEX, BER. TF 
算 模 型 本 文 不 作 有 共 体 描述 ， 详 情 见 参考 文献 心 201。 
Ek qg, 同时 受到 大 小 相同 ， 方 回 相 反 的 斥 力 ， 每 个 球 
体 所 受 的 总 斥 力 为 其 近邻 列表 内 其 他 球体 对 其 所 施 
FIREM. 
3.3 引力 及 弯曲 力 

在 便 核 系统 的 调整 过 程 中 ， 为 保证 纤维 构 型 不 
因 外 为 作用 下 过 上 度 扰动 ， 并 避免 过 度 弯 折 ， 在 同一 
纤维 球 链 的 球体 之 间 施 加 拉 压 弹 敬 和 这 曲 弹 车 ， 如 
图 7 所 示 。 


图 7 Jm ELS X B 


p Fa SE TR pfe De P en HH Fc, Hc 
弹性 常数 为 K,。 同 时 为 了 保证 球 链 的 紧凑 性 , 令 球 


q; 不仅 受到 相 邻 前 后 两 球 9 M qa WRS, 35 
受到 间隔 前 后 两 球 q,, qu, 的 弹簧 力 (不 包括 
纤维 两 端的 球体 ) 由 于 球 链 在 运动 过 程 中 前 一 个 球 
9 1 的 弹 赞 恢复 力 对 下 一 个 球 9 产生 作用 使 其 发 生 
位 移 ， 而 球 9 位置 的 改变 又 对 后 面 的 球 q, 产生 弹 
慢 恢 复 力 ， 如 此 循环 往复 下 ， 该 过 程 稳定 性 降低 且 
将 无 法 收敛 。 为 了 稳定 运动 ， 我 们 引入 了 一 个 摩擦 
因子 fuos ， 其 作用 类 似 于 摩擦 力 ,， 用 以 限制 纤维 
的 拉 压 弹簧 恢复 力 并 稳定 球 链 的 运动 ， 即 当 弹 得 因 
伸 长 或 压缩 长 度 变 化 量 在 0 n UAN, PERE 
恢复 力 ; 当 弹 区 长度 变化 量 逐 渐 增 大 到 O， m s 
筑 恢 复 力 的 强度 达到 正常 水 平 ， 如 图 8 所 示 。 其 数 
学 表达 式 如 下 : 


: «Uu (24) 
fs as (ô) = 1/ 2 e 2cos[(5| Ô n) / (ON Ox )] , Ob < 0 < Ô nax 
l ? Ô nax < 0 


摩擦 因子 


Ô vin Ó max 
长 度 或 角度 的 变化 量 


图 8 弹 千 摩 探 因 了 
定义 弹 昔 力 的 大 小 与 弹 营 长度 的 变化 星 呈 线性 
KR, FA r HIERA q; 受 球 链 上 下 一 个 球体 9 
的 拉 压 弹 千 恢复 力 定 义 为 : 


r/2—1q; —d;, 
F, 一 Kifs adna Ce ed| 


)(r/2—lg; =q) qi —i.i (25) 
la; -qin | 
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纤维 参数 进行 定义 0520。 对 于 本 文 的 模型 ， 式 中 摩 
探 因 子 中 的 0, 和 OO， 无 特殊 要 求 ， 选 取 条 件 较为 
宽松 ， 一 般 按 经 验 进 行 选取 ， 只 需 能 够 实现 算法 稳 
定 运 行 即 可 。 为 实现 本 文 的 模型 ， 其 分 别 取 为 0.01 
和 0.02. 除 此 之 外 , 球 9; 同时 还 受到 球体 9 5» di 
及 9 的 弹 和 恢 复 力 。 以 此 类 推 ， 可 以 得 到 整个 球 
链 系 统 所 有 施加 斥 力 后 的 球体 受到 的 拉 压 弹 敬 恢复 
Jj. 

在 球 链 保持 连续 性 的 同时 ， 需 要 保证 其 结构 的 
光滑 性 , 因此 在 球 与 球 间 施 加 角度 力 , BITS Hi SRL, 
RE ERO K, ， 同 样 根据 实际 纤维 参数 进行 定 
义 0520。 与 拉 压 弹 党 力 相 同 ， 每 个 球体 均 受 到 其 相 
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邻 四 个 球 的 角度 力 〈 不 包括 两 疹 球 )。 
对 球 9 受到 球 q; AER q, 的 角度 力 进 行 分 
析 ， 假 设 三 个 球体 初始 时 的 夹 角 为 元， 变化 后 夹 角 
Nyo Ika SIER q; ，、 球 9 ERREEN P, E 
义 角度 力 方向 为 9 p ， 大 小 与 球体 d 偏离 初始 坐标 
的 距离 成 线性 关系 ， 由 此 可 见 弯 曲 弹 筑 恢 复 力 即 为 
避免 受到 斥 力 的 重 二 纤维 球体 在 此 局 部 区 域 过 度 变 
形 ， 使 其 共有 回 直 纤维 形式 变形 的 趋势 。 同 时 考虑 
摩 探 因子 的 影响 ， 则 球 g 受到 球 q. MER q AES 
曲 弹 得 恢复 力 如 下 : 
F,, =Ksf; s CO-a) (26) 


ER q; KER q ,也 同时 受 球 dg; 的 角度 力作 用 ， 
方 回 与 相反， 大 小 为 其 一 半 。 同 理 , 球 g, 也 受 
FER qg KE qi WAE. 5 AME DCR 7J EET 
KAFERE d Hs REC 7B I FELI 
型 ， 摩 控 因 了 于 5 ， 与 O， 分 别 取 0.0028 及 0.0056. 
3.4 硬 核 系 统 调 整 的 实施 过 程 

上 述 算 法 首先 对 各 纤维 球体 遍历 进行 重 县 判 
汤 ， 己 经 明确 三 维 近 邻 列表 可 以 大 幅 加 速 这 一 判断 
过 程 ， 在 此 过 程 中 ， 某 一 纤维 球体 gq 可 能 同时 与 其 
近邻 列表 内 的 多 个 其 他 纤维 的 球体 同时 发 生 重 鞋 ， 
便于 算法 实施 , 选取 与 此 纤维 球体 重 琶 度 Al 最 大 的 
一 个 纤维 球体 q, 进行 记录 , 和 暂时 忽略 其 他 重 登 球体 
的 影响 。 随 后 ， 依 次 对 所 记录 的 球体 施加 斥 力 ， 再 
施加 拉 压 弹 车 恢复 力 及 过 曲 弹 党 恢 复 力 ， 且 两 种 弹 
早 力 作用 只 发 生 在 受 斤 力 的 球体 g, 自身 的 单元 组 
内 〔 即 同一 根 纤 维 上 与 球体 gq, 相 邻 的 四 个 球体 
Qio» Gas 4 及 9;,，)， 对 早 元 组 外 球体 的 影响 
不 再 做 考虑 ,以 倪 算 法 将 无 休止 地 进行 , 无 法 收敛 。 

所 记录 的 重 登 球体 g 所 受 总 力 即 为 其 近邻 列 
表 内 与 其 具有 最 大 重合 度 的 男 一 球体 对 其 施加 的 斥 
力 以 及 其 自身 单元 组 内 相 邻 球体 对 其 施加 的 弹 赞 力 
及 角度 力 的 矢量 和 。 在 对 记录 内 的 全 体重 有 登 球体 进 
行 上 述 操作 后 ， 再 次 进行 明 历 重 肢 判断 ， 重 复 以 上 
E. 

由 于 该 硬 核 调整 过 程 每 次 仅 针对 记录 内 的 一 对 
最 大 重 登 度 的 纤维 球体 及 其 目 号 单元 组 内 的 几 个 纤 
维 球体 ， 因 而 计算 速度 较 快 ， 也 能 保证 在 纤维 的 该 
局 部 区 域 具 有 物理 结构 的 真实 性 ， 同 时 整个 人 硬 核 调 
整 过 程 会 使 发 生 重 登 处 的 局 部 纤维 段 趋向 于 平 直 排 
列 ， 因 此 在 保证 较 高 的 纤维 体积 分 数 的 同时 ， 上 述 
过 程 也 容易 达到 收敛 。 

综 上 ， 通 过 对 外 部 周期 性 纤维 球 链 和 内 部 纤维 


球 链 分 别 记录 ， 施 加 上 述 作 用 力 ， 生 成 具有 周期 性 
几何 条 件 的 各 球 链 坐 标 ， 进 而 创建 包含 纤维 波纹 度 
特征 的 UD-FRP 高 保 真 RVE 模型 。 


4 RVE 模型 创建 及 探究 


4.1 三 维 高 保 真 RVE 模型 构建 

为 将 纤维 的 数字 模型 〈 球 链 上 的 球 心 坐标 及 各 
球 链 的 半径 ) 快速 转换 为 实体 模型 ， 本 文 基 于 
Microsoft Visual Studio 载体 进行 框架 搭建 〈 应 用 语 
AN CH), 并 适 配 建 模 软 件 CATIA V5 开发 集成 。 
通过 读 取 纤维 的 数字 模型 ， 对 纤维 进行 快速 的 实体 
建 模 ， 最 终生 成 具有 不 同 纤维 波纹 度 特征 的 
UD-FRP 三 维 实体 模型 。 

图 9 所 示 为 所 创建 的 不 同 纤 维 波纹 度 
( k, =200,400 及 直 纤 维 ) 下 ， 纤 维 体积 分 数 
58.9% ， 纤 维 半 径 45 um, RVE 尺寸 为 
108x108x108,um 的 RVE 几何 模型 示意 图 。 


k, = 200 


k, = 400 


直 纤 维 


9 不 同 纤维 波纹 度 的 RVE 模型 


4.2 纤维 波纹 度 对 横向 剪 切 性 能 的 影响 

上 文 创建 了 具有 纤维 波纹 度 特 征 的 RVE 几何 
模型 ， 选 用 T300/5284 复合 材料 建立 UD-FRP $f] 
剪 切 的 微观 力学 计算 模型 。T300 纤维 作为 FRP 广 
泛 使 用 的 一 种 增强 相 材 料 ， 领 域内 已 有 标准 性 能 参 
考 数值 20， 如 表 1 所 示 。 


机 械 工程 d 


表 1 T300 纤维 性 能 


参数 类 型 参数 值 
Ej/GPa 228 
E;-E;/GPa 17.2 
V12=V13 0.2 
V23 0.5 
G12=G13/GPa 27.0 
G23/GPa 5.73 


对 于 5284 树脂 基体 , 已 有 相关 研究 表明 其 对 毅 
KEJUR, KEH Drucker-Prager x78 $5835 
其 塑性 行为 ， 其 屈服 函数 如 下 : 

F =t— ptanf —d =0 (27) 
其 中 ，p 为 前 水 压力 ， NEÆERHM, d 为 树脂 基体 
的 内 聚 力 ，1 为 屈服 可 应 力 ， 均 可 由 下 式 得 出 : 


p=-trace(o)/3 (28) 
A (29) 
3- sino 
d ne (30) 
k 3 ! 
ERP 2 259019 ,为 三 轴 拉 伸 届 服 应 力 与 三 轴 压 
3+sin o 
缩 屈服 应 力 之 比 ，g 为 M-C 模型 的 摩擦 
sin o = 22 (31) 
O, +0, 


以 上 Go X o, 分 别 代表 材料 抗 拉 强度 及 抗 压强 度 ， 
届 服 前 应 力 f 则 可 由 下 式 表 示 : 
] 1 l.r 
- Sql —---2)-—» (32) 
(7 3dicc-ü- CY] 


其 中 ，g 代表 Mises 等 效应 力 ，r TUXETRINL7JUK 8 
的 第 三 不 变量 。 

5284 树脂 在 承受 里 轴 拉 伸 载 向 时 , 在 相对 较 小 
的 应 变 下 即 会 有 友 生 脆性 断 委 ， 在 单 轴 压 绾 或 筋 切 载 
何 下 ， 则 上 其 有 明显 的 届 服 及 塑性 变形 。 为 准确 表示 
树脂 材料 在 不 同 条 件 下 的 损伤 行为 ， 采用 Ductile 
Damage 定义 树脂 基体 的 损伤 起 始 及 演化 ， 基 体 刚 
度 随 损伤 起 始 条 件 的 不 同 而 降低 ， 性 能 标定 如 表 2 
Bran o 


表 2 5284 树脂 基体 性 能 


参数 类 型 参数 值 

弹性 模 量 E/GPa 1.8 

拉 伸 屈服 强度 co/MPa 17.34 
压缩 届 服 强度 c, MPa 82.75 
抗 拉 强度 O, /MPa 58.78 
抗 压 强度 ca./MPa 160.28 
Drucker-Prager 摩擦 角 / B 47.63* 
流动 应 力 比 /K 0.778 


断裂 应 变 /6Ev 0.03078 
应 变 能 释放 率 /C 0.2J/m? 


针对 含有 不 同 纤维 波纹 度 k, 的 RVE 力学 模型 
展开 横 回 剪 切 的 有 限 元 计算 。RVE 的 宏观 力学 性 能 
计算 对 RVE 尺寸 及 网 格 大 小 均 具 有 一 定 的 敏感 性 ， 
因此 有 必要 进行 相关 收 钱 性 分 析 , 据 现 有 研究 结果 ， 
针对 UD-FRP 的 RVE, 其 三 维 尺 寸 均 为 纤维 半径 的 
24 f CHI RVE 尺寸 为 24rx24rx24r ) 时 ， 网 格 尺 
寸 为 纤维 半径 的 va 倍 《〈 即 网 格 尺 寸 为 r14) 时 ， 
有 限 元 分 析 的 计算 结果 有 效 ， 可 被 用 于 宏观 力学 性 
能 分 析 07。 本 例 中 采用 的 T300 纤维 半径 为 4.5 um ， 
因此 利用 上 文 所 述 算法 ， 创 建 k, =200， 三 边 尺寸 
为 24:r， 即 108x108x108um 的 RVE 模型 ， 由 于 
纤维 与 基体 各 自 的 几何 形状 较为 复杂 ， 采 用 C3D4 
三 维 线性 四 面体 进行 网 格 划 分 ， 网 格 尺 寸 为 
r/4-1.125um, WE] 10 所 示 。 


Face 2 


Face 4 b. 
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基于 RVE 建 模 的 一 个 重要 方面 是 RVE 的 尺寸 
大 小 和 施加 的 边界 条 件 ， 应 用 的 边界 条 件 将 会 影响 
整体 的 机 械 性 能 P2。 此外, 当 有 限 元 模型 足够 大 时 ， 
可 以 忽略 边缘 和 面 效 应 多。 因此 本 文采 用 标准 边界 
条 件 , 将 直接 约束 应 用 于 RVE 的 边界 , 即 运动 统一 
边界 条 件 (KUBC)， 具 体 设 置 如 下 : 

面 2 沿 全 厚度 〈y 方 同 ) SAHE, By 
位 移 y FW) 应 用 于 面 1, 在 面 3 和 面 4 之 间 应 用 
连接 约束 (Tie Constraints)， 面 S5 和 6 可 以 目 由 变 
JE. 

有 限 元 计算 所 得 RVE Bf [93] 898 77] Iz 73 MA HH 
线 如 图 11 所 示 ， 结 果 初 步 表 明 ， 不 同 的 后 值 下 ， 

曲线 的 和 斜率 未 见 显 闭 差 别 ， 即 纤维 波纹 上 度 特征 对 
UD-FRP 的 横 同 雯 切 模 量 没 有 明显 的 影响 ， 然 而 纤 
维 波纹 度 特征 对 其 横 回 驴 切 的 强度 有 着 一 定 的 强化 


We BHSE: 


EH, ARRARRK BIK HER), MJN 
变 曲 线 的 峰值 越 高 ， 即 UD-FRP 的 横向 剪 切 强度 越 
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0 
0.000 0.005 0.010 0.015 0.020 0.025 
fs E BIYA vy. 
11 不 同 纤维 波纹 度 下 UD-FRP 横向 剪 切 应 力 应 变 曲 线 图 


这 种 影响 规律 初步 竹 夸 是 由 于 纤维 的 雪 曲 使 得 
纤维 与 基体 在 任意 一 个 忽 纹 扩展 面 内 的 接触 面积 增 
加 ， 因 而 纤维 为 基体 承担 了 更 多 应 力 ， 同 时 也 更 加 
容易 阻碍 基体 神 纹 的 扩展 ， 因 此 需要 更 大 的 应 力 太 
应 变 才 会 导致 前 切 损 伤 破 坏 的 友 生 ，UD-FRP 的 模 
[3] 95 1) o8 JE fa PA en o 


5 th 论 


本 文 针 对 单 癌 纤维 增强 复合 材料 IUD-FRP)， 基 
于 纤维 随机 游 走 算法 ， 结 合 三 维 近 邻 列表 ， 在 提高 
算法 速率 及 准确 性 的 同时 ， 保 证 了 代表 性 体积 单元 
CRVE) 的 几何 周期 性 不 被 破坏 ， 成 功 构建 出 可 控 
纤维 波纹 度 的 高 保 真 几何 结构 。 随 后 针对 创建 的 模 
型 分 别 展 开 了 横 癌 前 切 计 算 ， 得 到 以 下 结论 : 

(1) 通过 多 变量 妆 米 斯 - 费 歇 尔 分 布控 制 纤 维 
的 波纹 度 特征 ， 在 调控 波纹 度 大 小 的 同时 ， 通 过 该 
分 布 自身 的 特性 保证 了 纤维 的 物理 真实 性 ， 避 人 免 了 
纤维 在 局 部 产生 弯 折 过 度 的 问题 ， 更 加 贴 合 实际 的 
纤维 结构 。 

(2) 改变 算法 中 的 可 信和 度 参 数 太 及 态 即 实现 
RVE 模型 纤维 波纹 度 的 更 改 , ,设置 越 小 则 所 建 纤 
维 的 波纹 度 越 高 , 最 终 通过 CATIA 的 二 次 开发 技术 
进行 实体 模型 的 快速 创建 。 

(3) 通过 对 不 同 纤 维 波纹 度 的 RVE 展开 模 回 
前 切 强 度 的 有 限 元 计算 ， 对 比 计算 结果 可 以 初步 看 
出 ， 这 种 波纹 度 特征 会 对 UD-FRP 的 横向 剪 切 强度 
造成 一 定 影响 ， 纤 维 波纹 度 越 大 ，UD-FRP jf 
剪 切 强度 越 大 。 
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